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RANGE-ENERGY  MEASUREMENTS  FOR  742-Mev  PROTONS 

John  E.    Ward 

Lawrence  Radiation  Laboratory- 
University  of  California 
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May  4,    1959 


ABSTRACT 

This  is  a  proposaj.  for  an  experiment  to  measure  the  relative 
stopping  powers  of  Be,    Al,    Fe,    Cu,    Ag,    W,    Pb,    and  U  by  magnetically 
analyzing  the  degraded  proton  energies  of  the  742-Mev  proton  beam 
of  the  Berkeley  cyclotron  after  passage  through  the  absorbing  elements. 
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I.      INTRODUCTION 

If  a  charged  particle  passes  through  a  substance  with  sufficient 
energy  that  the  valence-shell  effects--i„  e„  ,    atomic  bindmg--can  be 
neglected,    the  average  rate  of  energy  loss  due  to  ionization  only  is 
expressed  by  the  well-known  Bethe  equation, 


dE         4.e4z2        NZ 


dx  2 

mv 


in     2mv_2    _  in(l-(32)-p2 


(1) 


where     e     and    m     refer  to  the  electronic  charge  and  mass,    ez  is  the 
charge  of  the  incident  particle,    NZ  is  the  number  of  electrons  per 
unit  volume  of  the  stopping  element,    and    I    is  the  "mean  excitation 
potential"  of  the  atoms  in  the  stopping  element. 

The  development  of  this  formula  makes  use  of  the  Born  approxi- 
mation,   which  requires  the  orbital  velocity  of  the  electrons  to  be  small 

compared  with  the  velocity  of  the  incident  particle.     A  proton  with  a 

2 

kinetic  energy  greater  than  2  Mev  satisfies  this  condition. 

It  can  be  seen  that  the  stopping  power  depends  only  on  the 
velocity  and  charge  of  the  incident  particle  and  not  its  mass,    and  the 
constant  of  primary  importance  as  far  as  the  stopping  substance  is 
concerned  is  its  mean  excitation  potential.      The  mean  excitation 
potential  is  defined  as  the  product  of  all  excitation  levels  in  the  sub- 
stance,      fico.,  ,    each  weighted  exponentially  by  its  respective  oscillator 
strength,    f.,  :    ° 

fik 
I  =  KlTcoik  ,  (2) 
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where 


k  lk 


This  formula  would  offer  a  means  of  calculation  of  I,  but,  except  for 
the  lightest  of  elements,  a  very  impractical  one.  I  is  usually  taken 
as  an  experimentally  adjustable  constant. 

Bloch  has  shown,    using  the  statistical  model  of  the  atom,   that 
in  the  Thomas -Fermi  approximation  one  has 

I  =  KZ,  (3) 

3 

where    K    is  an  adjustable  constant.  "     Using  the  more  accurate  Thomas 

4 
Fermi-Dirac  model,    Jensen  found 

IQ  =  K0(l+k0Z-2/3)  Z.  (4) 

The  index  0  refers  to  isolated  atoms.      This  function  is  plotted  in 
Fig.    1  as  the  solid  line,, 
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Fig.    1.      Experimental  mean  excitation  potentials  and  the 
statistical  atomic  model.     O:  I/Z  of  bound  atoms  meas 
ured  with  340-Mev  protons  and  evaluated  relative  to 
I.  .  =  165  ev;   #:  Iq/Z  of  isolated  atoms.     A:   i/Z  meas- 
ured with  10-20-Mev  protons  and  corrected  for  inner 
shells.      The  corrections  apparently  underestimate  the 
stopping  -power  deficiencies  of  inner  shells  at  low 
energies. 


IL     EXPERIMENTAL 

Since  at  high  energies  the  corrections,  applied  to  the  Bethe 
formula  for  valence-shell  effects  are  negligible,    the  742-Mev  protons 
from  the  Berkeley  cyclotron  seem  to  be  a  natural  choice  for  the  incident 
particles.     In  addition,   this  energy  range  would  serve  as  a  check  on  the 
340-Mev  Bakker  and  Segre  data. 

Figure  2  is  a  reproduction  of  a  slide  used  by  Werner  Brandt  at 
the  Gatlinburg  Meeting  on  Penetration  (September,    1958),    and  compares 
the  effective  mean  excitation  potentials   (without  corrections),    meas- 
ured by  Burkig  and  MacKenzie  at  UCLA  in  1957,      with  the  fluctuations 
of  the  electron  densities  in  the  periodic  system,,      One  can  notice  slight 
but  significant  trends  in  the    I    values,    which  are  parallel  to  the  electron 
density  wherever  the  accuracy  in  I    is  high  enough.      The  second  max- 
imum should  have  been  shown  clearly  by  Mo (42)  and  Rh(45),    but  un- 
fortunately these  were  the  poorest  foils  and  hence  gave  the  least  accur- 
ate data  in  the  series.      Lead  showed  an  "anomalously"  high  value  and 
needs  an  independent  check.      Thorium,    on  the  other  hand,    confirmed 
the  expected  trend. 

On  the  basis  of  this  diagram  and  on  the  advice  of  Dr.    Brandt 

the  following  metals  were  chosen  as  target  elements:     Be (4),    Al(13), 

Fe(26),    Cu(29),    Ag(47),    W(74),    Au(79),    Pb(82)s    and  U(92).     It  was 

decided  to  use  enough  absorber  to  degrade  the  protons  approximately 

15%  in  momentum  or  from  742  Mev  to  570  Mev  in  kinetic  energy.      In 

each  case  an  equivalent  amount  of  absorber  was  to  be  used  so  that  the 

results  could  be  compared  on  a  relative  basis  with  aluminum  as  the 

standard.      (See  Table  I)    Adequate  information  exists  on  Cu,    Al,    Ag, 

Pb,    and  Be  to  determine  the  amount  of  absorber  required.      Specifically, 

7 
the  tables  and  graphs  m  UCRL-2426  (rev.  )     were  used.      To  determine 

the  equivalent  thickness  of  Fe,    W,    Au,    and  U  a  numerical  integration 
of  Bethe'  s  formula  was  performed  using  values  of    I    previously  de- 
termined by  Bakker  and  Segre  (see  Appendix  A).      Clearly  these  values 
are  not  precise.      To  provide  for  a  means  of  interpolation  in  order  to 
use  the  precise  equivalent  amounts  of  absorbers  for  the  relative 
comparison,    an  additional  length  of  aluminum  absorber  is  to  be  added 
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Fig.    Z.      Mean  excitation  potential  compared  with  electron 
density.      (From  Werner  Brandt,    Gatlinburg  Meeting 
on  Penetration,    Sept.    1958.  ) 
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to  each  absorber  in  order  to  provide  an  additional  point.      This  absorber 
"shim,  "  9  g/cm     of  aluminum,    is  to  be  placed  first  on  the  "upstream" 
and  then  "downstream"  side  of  the  various  absorbers  and  the  results 
averaged  to  give  the  effect  of  adding  this  additional  absorber  in  the 
middle.      The  equivalent  values  of  the  aluminum  shim  are  listed  in 
Table  II. 

Table  I 


Equivalent  absorber 

thicknesses 

Absorber 

Density 
(g/cm    ) 
1.85 

Length 

(g/cm    ) 
85 

(cm) 

45.94 

(in.) 

Be(4) 

18.09 

Al(13) 

2.70 

89 

32.98 

12.98 

Fe(26 

7.8  5 

91 

11.59 

4.56 

Cu(29) 

8.93 

101 

11.30 

4.45 

Ag(47) 

10.49 

110 

10.49 

4.13 

W(74) 

18.60 

115 

6.20 

2.44 

Au(7  9) 

18.88 

118 

6.24 

2.46 

Pb(82) 

11.34 

133 

11.73 

4.62 

U(92) 

18.70 

125 

6.68 

2.63 

Table  II 

Equivalent  value 

of  Al  shim 

Absorber 

3- 

Equiv.    value 

Thic 

kness   (g/cm    ) 

Be(14) 

1.024 

9.216 

Al(13) 

1.000 

9.000 

Fe(26) 

0.906 

8.154 

Cu(29) 

0.87  5 

7.875 

Ag(47) 

0.789 

7.101 

W(74) 

0.680 

6.120 

Au(79) 

0.676 

6.084 

Pb(82) 

0.660 

5.940 

U(92) 

0.630 

5.670 

a 

mass  sto 

pping  power 
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The  experimental  setup  is  shown  in  Fig.    3.      The  full  energy  of 
the  184-inch  cyclotron  is  used,    742  Mev.      The  proton  beam  passes 
through  a  pair  of  jaws  reducing  it  to  l/l6  in.    in  the  horizontal  dimen- 
sion.     The  beam  then  passes  through  focusing  and  bending  magnets 
and  comes  to  a  focal  point  1  foot  within  the  physics  cave.      This  point 
is  the  physical  center  of  the  absorbers.      The  beam  will  then  pass 
through  a  focusing  magnet  (three-element  quadrupole),    and  a  bending 
magnet  in  which  the  beam  will  be  momentum -analyzed,    and  then  fall 
upon  x-ray  film  placed  at  the  focal  point  of  the  quadrupole. 

The  procedure  will  be  to  obtain  an  x-ray  photograph  of  the  full 
beam  and  then  a  photo  of  the  beam  with  the  absorber  in  place  for  each 
of  the  elements  listed.     In  addition,    photos  will  be  taken  with  the 
aluminum  shim  as  previously  explained. 

The  density  of  the  developed  x-ray  film  will  be  proportional  to 
the  number  of  protons  striking  it.      By  reading  the  film  with  a  den- 
sitometer one  can  obtain  a  plot  of  the  numbers  of  protons  versus  their 
positions  on  the  film.     In  oxder  to  facilitate  density  measurements  a 
duplicate  of  each  of  the  exposures  mentioned  above  will  be  made  at 
precisely  one -half  the  time  for  the  original.     Then,    for  example,   the 
maximum  intensity  on  the  second  exposure  will  be  exactly  one -half 
that  of  the  maximum  intensity  on  the  first,    and  can  be  used  as  the 
reference  for  finding  the  points  of  half-maximum  intensity  on  the  first 
film.      This  procedure  can  then  be  extended  ad  infinitum  throughout 
the  intensity  spectrum.      The  plot  of  numbers  of  protons  versus  film 
positions  is  expected  to  look  like  Fig.    4.     In  addition,    by  the  use  of 
wire  orbits,    a  plot  of  energy  versus  position  on  film  can  be  made.     This 
plot  should  look  something  like  Fig.    5.      The  two  graphs  can  then  be 
combined  to  form  a  plot  of  number  of  protons  versus  energy.     An  ex- 
ample of  this  type  of  plot  is  shown  in  Fig.    6.      In  general  it  has  the  same 
shape  as  Fig.    4.      This  plot  (Fig.    6)  will  be  the  tool  for  describing  the 
experimental  results. 
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Fig.    3.      Experimental  setup,    showing  placement  of  magnets. 
Not  to  scale. 
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Fig.    4.      Expected  general  contour  of  plot  of  the  number  of 
protons,    i.  e.  ,    the  film  density,    versus  position  on  an 
x-ray  film  after  undergoing  magnetic  analysis.      Two 
curves  are  shown:     that  of  the  full-energy  proton  beam 
of  the   184-inch  cyclotron,    74Z  Mev;   and  that  of  the 
same  beam  after  approximately  15%  momentum 
degradation. 
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Fig.    5.      Expected  general  contour  of  the  results  of  a  wire 
orbit  analysis  of  an  analyzing  magnet. 
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Fig.    6.      Expected  results  of  the  combination  of  data  ob- 
tained from  figures  4  and  5  showing  the  energy  spectrum 
of  the  full-energy  proton  beam,    742  Mev,    and  the  same 
beam  after  approximately  15%  momentum  degradation. 
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III.     DATA  ANALYSIS 

For  determination  of  the  mean  energy  loss  through  each 
absorber  the  center  of  gravity  of  the  two  curves  of  Fig.    6  must  be 
found.      This  can  be  done  easily  for  the  first  curve,    representing  the 
full  beam  of  protons  undegraded.     However,    because  of  close  collisions 
of  protons  and  electrons  the  second  curve  has,    theoretically,    an 
infinite  tail.     Obtaining  this  infinite  tail  is  not  practicable  experimentally, 
so  that  it  must  be  supplied  by  a  theoretical  calculation.      The  observed 
ranges  of  individual  particles  from  any  monoenergetic  source  will  show 
a  substantially  normal  distribution  about  the  mean  range.      Because  the 
very  hard  collisions  are  few  in  number,    the  actual  distribution  is  some- 
what asymmetric,    with  a  long  tail  in  the  direction  of  the  short  ranges, 

Q 

and  a  mean  range  which  is  slightly  shorter  than  the  modal  range.        For 
the  approximation  of  this  tail,    the  collisions  suffered  by  the  protons 
with  the  electrons  of  the  absorbers  will  be  divided  into  two  classes, 
"near"  and  "far„  "    The  so-called  near  collisions  are  those  responsible 
for  the  occasional  large  energy  losses  and  the  far  collisions  are  those 
contributing  only  to  the  mean  energy  loss,    but  introducing  negligible 
energy-loss  fluctuation.     In  other  words,    each  proton  will  suffer 
approximately  the  same  minimum  energy  loss  and  may  in  addition 
suffer  a  near  collision  (or  collisions)  resulting  in  a  much  greater  energy 
loss.     Each  absorber  will  then  be  divided  into  a  number  of  equal  seg- 
ments in  which  the  probability  of  a  near  collision  is  known,    the  energy 
loss  will  be  computed  for  this  one  segment,    and  the  resulting  spectrum 
will  be  "folded"  into  itself  as  many  times  as  there  are  segments.      It 
is  hoped  that  this  spectrum  will  fit  the  experimental  curve  and  provide 
the  needed  tail.      This  procedure  seems  valid  in  that  the  so-called   near 
collisions    are  those  in  which  the  binding  energies  of  the  electrons  are 
insignificant,    and  this  is  the  assumption  made  originally  by  neglecting 
the  corrections  to  Bethe'  s  formula. 

By  this  means  of  completing  the  graph  for  the  degraded  protons, 
the  center  of  gravity  of  the  curve  can  be  found  and  the  mean  energy  loss 
determined.      This  procedure  would  be  done  for  each  of  the  elements. 
By  interpolation  between  the  two  curves,    obtained  by  use  of  the 
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aluminum  shim,    for  each  of  the  elements  the  thickness  of  absorber 
equivalent  to  aluminum  will  be  determined.     The  relative  mass  stopping 

power  of  each  element  will  be  calculated  by  dividing  the  amount  of 

2  ? 

aluminum  (in  g/cm    )  by  the  amount  of  the  element  (in  g/cm    )  needed 

to  produce  the  same  mean  energy  loss. 

In  addition,    the  stopping  power  per  electron  is,    from  Bethe'  s 

formula, 

dE 


4    2 
dX  4ue    z 


NZ  2 

mv 


Cn^L.    -in(l-p2)-|32 


The  relative  stopping  power  per  electron  with  respect  to  aluminum  is 


q  = 


-  inl  +  [in2mv2  -  in(l -(3  2)-(3Z] 

-  *nIA     +  [in2mv2-in(l-p2)-p2] 


This  quantity  can  be  easily  calculated  by  multiplying  the  relative  mass 

stopping  power  by 

A  zAi 

X     -~—     . 


AAi 


Thus     the  means  for  calculating    I    is  apparent.     In  this  calculation  the 

9 
value  of  I  =  150  ev  for  aluminum  is  used. 

From  the  values  of    I    calculated,    a  curve  of  i/Z  versus     Z    will 
be  plotted.      This  curve  will  in  general  look  like  Fig.    1.      The  value  of 
the    Bloch  "constant"  can  then  be  compared  with  previously  obtained 
values. 

In  order  to  determine  the  length  of  absorber  traversed  by  the 
protons  a  correction  for  multiple  scattering  must  be  made.      In  the 
usual  approximate  correction,    the  observed  range  is  expressed  as  a 
projection  of  the  mean  range: 

R   .        =     Y   i.  cos0.    -       V    !,  (1   -  I  QZ) 
obs  i—     \  i  {_      i  2     i 


=  r  -  \     y  i.    ez 

mean        Z        Z_     i        i 
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where        i.     is  the  distance  between  the  ith    and  the  (i-l)st  small- 
1  — 

angle  collision   and  6.  is  the  direction  with  respect  to  the  beam  axis 
after  the  ith  collision,,      Therefore  we  have 

R 


obs 

R  -  R   ,      =    AR  =  i- 

mean        obs  2 


<  es2'x»> 


dx 


where   \0~    (x)  \     is  the  mean  square  deviation  in  angle  from  the 
normal  direction  and    x    is  the  projected  distance  in  the  absorber. 
Also, 


0  =      J  6c2(x)) 

rms  \ S        / 


21  Mev 
(3cp 


>/T 


where    t    is  the  absorber  thickness  measured  in  radiation  lengths. 
Using  the  value 
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(e2(x)> 


R 


0 


8  TT  Ne4Z2G(pv)~2dx> 


x 


Bichels    Mozley,    and  Aron  found 


13 


AR 
R 


=  (a -R.kJA 


mean       obs  '      mean 


Gm.c    Z 
(M0c2)2 


xEQ         ^2  rEQ      _1(32 

'o  b7z  dKk     pWz)    dE 


£ 


b/z 


dE 


Where  G=      2i 


-1/3 
n  181   Z      '  \    B  i 


s  the  stopping  number  of  the  atom, 


B  =  Z 


in 


2mv 
I(l-(32) 


-  P' 


and  mn  is  the  rest  mass  of  the  electron.      The  ratio  AR/R  is  a  slowly 
varying  function  of  the  energy,    amounting  to  several  percent  for  high  Z. 
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The  corrections  found  by  these  experimenters  appear  in  Table  III.     In 
each  case  it  is  clear  that  as  Z  increases  the  correction  for  multiple 
scattering  increases,    and  as     E    increases  this  correction  decreases. 

Table  III 


Multiple -scattering  corrections 


Element  E 


Range    , 
6   ob 


a 


(Mev)  Path  length 

Be  17.340  0.9990 

9.588  0.9988 

Al  17.836  0.9962 

14.971  0.9960 

11.820  0.9956 

6.150  0.9943 

Cu  17.893  0.9900 

9.938  0.9872 

Ag  17.923  0.9818 

10.022  0.9759 

Au  17.549  0.9626 

9.698  0.9475 


correction  for  multiple  scattering 
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APPENDIX 

A.     Determination  of  Absorber  Thicknesses 
1.      Fe,    W,    Au,    and  U 

Since  no  usuable  data  exist  on  these  elements,    a  numerical 
integration  of  Bethe'  s  equation  was  performed  in  order  to  determine 
the  absorber  thickness.      The  total  desired  energy  drop  was  divided 
into  twelve  equal  steps  such  that  the  change  in    (3     was  constant, 
0.004.      The  energy  increment  was  then  calculated  f  rom 


E  =  m0Y°    ' 


where 


Y  = 


\l  1  -f3z 


Therefore  we  have 


AE  =  mQc2(3  y2  A(3  =  3.752  [3  y3  (Mev). 

The  integration  on  Bethe'  s  equation  was  performed  by  de- 
termining 


rdE 

^dx   ; 


ZBpZ 
(3A 


^GP2y2 


G    J 


A(3, 


where 


B 


4    2 
4Tre    Z    Nq 

2 
mQc 


=     7.21x10 


-7 


g  cm 

T 
sec 


atoms 
g-atom 


G  =    h^L     x  106. 


The  values  determined  appear  in  Tables  IV  and  V.      Graphs  of  the  re 
sultant  variation  in  dE/dx  are  shown  in  Figs.    7  through  10. 
2.      Be,    Al,    Cu,    Ag,    and  Pb 

Sufficient  data  on  these  elements  exist  in  UCRL.-2426  (rev). 
3.      Summary 

The  absorber  thicknesses  to  be  used  were  given  in  Table  I. 
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Fig.    7.      Calculated  range-energy  curve  for  Fe. 
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Fig.    8.      Calculated  range-energy  curve  for  W. 
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Fig.    9.      Calculated  range -energy  curve  for  Au. 
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Fig.    10.      Calculated  range -energy  curve  for  U. 
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B.      Bending -Magnet  Analysis 


1.     Angle  of  Bending 

If  the  Hi  of  the  magnet  and  the  Hp  of  the  particle  are  known, 
the  following  is  the  derivation  for  the  angle  of  bending  (See  Fig.    11): 


6 


sin  j 

2p"       ' 

= 

(Hp)        ..    , 
r  particle 

P 

H 
magnet 

".    sin 

e 

2 

(Hi) 

mag, 

2<H<»part. 

(Hp) 


part. 


(Hi) 


mag. 


X     i; 


a  a 

if    sin  —    ^  then    0 


(Hi) 


mag. 


(Hp) 


part. 


This  approximate  formula  works  well  for      G    less  than  1  radian. 


Fig.      11.      Path  of  a  particle  in  a  magnetic  field. 
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2.     Vertical  Focusing  in  an  Analyzing  Magnet 


y 
X 


>\ 


z    «*-«*l    Ij    j) 

*-*/*V 


K 


L 


Fig.    12.     Side  and  top  views  of  typical  bending  magnet,    showing 
field  configurations  and  entrance  and  exit  angles. 


In  Fig.    12,    let  0     =  d     =    |     and  tanG^-  -J£ 


From  Maxwell'  s 


equations  we  have 


V  X  H     =  0S 


3H  9H 

y  x 

3x  — <Ty~ 


-  0, 


6H 


3H 
i 

9z 


=  0, 


9H 


x 


dz 


9H 

z 

"97 


(1) 


=  0; 
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v  •   H     =    0, 


3H  3H  9H 

-x  +      y  +  _^  =  o. 


3x  3  y  3  z 


F     =     (e/c    v   X  H) 

y  y 


=  e/c(v  H    -v  H   )  =  -e/c  v    H   . 
'  Z     X         x     z'  '  X      z 


3H 

Assuming      — ^ —         is  constant  leads  to 
0  3  v 


y 

3Hz  3Hz  3H 

a  Ay  =         .  y    =    -jr- ^     y. 

z         /  3y  '  3y        '  3z        ' 

0 


Therefore 


3H 

F     =  -  e/c  v    H     =  -  e/c  v      — *-*-    y. 
y  x     z  '         x       3  z       ' 


I  shall  let    H     =  0;  then 
x 

3H  3H  3H 

and  substituting  Eqs.    (2)  into  (1),    one  has 

3H  3H  3H  3H  3H  3H 
7                z    .  fl             z             _y             _y_  __z 

3x  3x         u'        3y  3z    »       3y  3z      * 

Using  Gaussian  units, 
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The  momentum  imparted  to  the  particle  in  the    y    direction  is 

Ap    ==    fFdt=[F     ^    =   .1  (Jjl   y    fjr    dz 

y    J         j   y  vz        c  J  vz       dz 


aH         \ 

=    -     tan   e.y  U-i      dz  J     =    -    tanG.yAH 
c  1  '  \  dz  /  c  1  7         y 


Defining  6  =  Ap   /p,     where    p  =  eHp/c,    one  has 


(e/ckane,  yAH  AH 

'                     1                 V  ^         V                   V 

6  =        /    /   vtt *-  =    tan  0.    ■*-         „  y 

(e/cjHrp  1     p          H 


In  the  pole  gap    H  =  H    ,    and  AH     is  the  increase  of  H     from  a  place 
of  zero  field  to  the  edge  of  the  pole  tips,    where  H     becomes  uniform 
and  equal  to  H.      Therefore 


and 


AH     =  H 

y 


6    =  tanG,     ^    .  (3) 

1     p 
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I  have  assumed  H     small  until  reaching  the  pole  tips,    so  that 
v     is  constant.      From  Fig,    13  we  have 


y  =  S0(3  =  S^. 


y 

Z     f^     X 


Fig,    13,      Vertical  bending  at  pole  tips  due  to  field  variations. 


Therefore 


or 


H       S.  SQ 


1         1  6  1_ 

S0  +  S.-      y  f     ' 


from  Eq.    (3), 


1_     _    ^  tan_#i 

f  y  p 

This  means  that     f  'is  independent  of  the  height  at  which  the  particle 
hits  the  entrance  of  the  magnet.      This  independence  contributes  the 
focusing  action  in  the  vertical  plane.      The  principal  planes  are 
assumed  to  be  the  ends  of  the  magnet.      In  terms  of  the  total  angle  of 
bending. 

1  tan  0/Z 


T 


(4) 
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3.      Change  in  Vertical  Focal  Length  with  a  Change  in  Momentum 


Given  f  =    ^   P    -       ,  sin0,   = 


tan   #j  1  2P 

we  have  n  . 

sinfl,  1 

tan  9,   - 


and 


or,    using 


si  1  -  sin20, 

1 


Therefore 

f  = 

P 
P 

f 

"J 
-J 

1    4p2 

v  sin   9, 

A 

PC 

eH    » 

1 

-         -      1 

L2 

since 

/4p4 

2 
-   P 

Af  =   f     \*4     "  2 1    A<" 


Ap  =  (c/e  B)  Ap  =  p  Ap/p, 


Af  =  ¥r   ( if¥r  -'I'p. 

e    H   f    leTL' 


A£  =    i!  (iE2.    .A    ^p_    _  (5) 
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4.      Dispersion  in  an  Analyzing  Magnet 

In  order  to  determine  the  dispersion  particles  of  various 
momenta  passing  through  an  analyzing  magnet,    it  is  necessary  to  find 
distance    d    in  Fig.    14.      Figure   14  shows  two  particle  tracks,    with 
point  (h,  k)  being  the  center  of  curvature  of  the  lesser  momenta 
particle.      From  the  geometry  we  have 

h  =  -  L/2  +  p1    sin(e/2), 
k  =  p<  cos(0/4 
The  equations  for  the  circles     s    and    sr    are 
(x-0)2  +  (y-pcos(<9/2))2  =  p2, 

(x+  \  -  p'  sin(e/2))2  +  (y-  p«  cos(<9/2))2  =  p<  2. 

To  find    d    it  is  then  necessary  to  find  the  intersection  of  the  circle 
s1    and  the  line  x  =  L/2.     Note  that  we  have 


sin(0/2)=  l^t  cos(e/2)=    ^  ^  '  L^4      . 

Therefore  the  equation  for  the  circle  s "    is 
[x+(p-p')sin(0/2)]2+  (y  -  p'  cos(<9/2))2  =  p»  2 

=  x2+2x(p-p'  )sin(0/2)+  (p-p»  )2  sin2(0/2)  +  y2  -  2yp«  coste/2)+(p'  )2cos2f0/2)=  p' 

=  x2+  2x(p~p'  )  sin(0/2)+  p2sin2(fl/2)-  2pp'  sin2(0/2)+y2  -  2yp<  cos(6>/2)=  0. 

Now  at  x  =  L./2  (i.  e.  ,   the  intersection  of  circle  s'    and  line  x  =  L./2), 
we  have 
2 
^-    +L(p-p'  )sin(0/2)+  p2sm2(0/2)- 2pp'  sin2(6»/2H  y2  -  2yp'  cos(G/$=  0. 

This  reduces  to 


y  =  p    cos 


(0/2)  -  n/P,2cos2(0/2)-   L2(l  -p'/p)         • 


33- 


(h,kv 


MU-  17452 


Fig.    14.      Dispersion  in  an  analyzing  magnet  of  charged 
particles  of  different  momenta. 
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Observing 


0         L 


sin2  = 


2p 


sm 


2p'      '        2'    "      2  2      ' 


we  see 


sin- 


.      0       60 
Sln   2   +  ~ 


L'  0  60  0  60 

-= — r    =    Sin^r  COS  tt—    +  COS  rr  sm-^- 

2p*  2  2  2  2 


0   t  60  0 

*    sin2  +  —    COS2 


Further, 

2p 

60             0        L1 
T   cos  2=^p"« 

2p- 

L 
2p 

L 


\li+d'/; 


2p'  2p 


Therefore 


■f    cos  I    60  =    A     N/l  +  d2/:L2      -  1 
L.  2  ps 


Since       p  =  L/2  siny 


P 
P1 


2     \ 


1  - 


2L 


-   1, 


we  have 


2  sin-=-    L 


y        cfl  P-P     l      P 

COS-r        60    =     i-— »-    +      •*- 


2L 


and 


60  =  2tan-| 


P-P    +   £ 


I    P*  P*       2L"J 
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Since    p  oc   p 


we  have 


PJ  P' 


and 


50  =  2tan- 


P-P'  +    P       _d_ 


2L 


where 


d  =  p'<    cos-^-    -     \f  p»      cos      -    -  L    (1  -  p/p?  ) 


Therefore 
60 


Pl      COSy 


■  ^4  (*£ 


+  2    p>  V     p1; 


.       2    0 
tan     ^ 


The  assumptions  made  were 


.     0  60  .0 

sin-=-  cos  -=—  =  sinj    , 


60  0        60  0 

Sin— j     COSy    =    ■=—       COSy        , 


/ 


A1  A1 

1   _L    d  -     H     4.       d 

,  +  _     1    +      

i/  2l/ 
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Wire  Orbits 


1.      The  Wire -Orbit  Method 


This  is  an  experimental  method  for  determining  the  path  of  a 
charged  particle  in  a.  magnetic  field.     It  is  based  on  the  analogy  be- 
tween the  shape  assumed  by  a  current-carrying  wire  in  a  magnetic 
field  and  the  path  of  a  moving  charged  particle  in  the  same  field. 

The  motion  of  a  charged  particle  in  a  magnetic  field  is  described 
by  the  equation 

2  ^  ^ 

=    ev     x   B 


mv 
~R~ 


with  the  solution 


BR  =  p/e. 


(1) 


The  forces  on  the  wire  in  Fig„    15  consist  of  the  tension,      T, 
and  the  force  due  to  the  magnetic  field, 

F  =  i  di     X   B  . 


T  * 


Fig.    15o      Forces  involved  in  the  wire -orbit  method- 


Equating  the  forces. 


2T  sin0  =    i  di  B, 


and  noting 


sin  B    -  — ^ —  as  di  -*-  0S 

XV 


one  finds 
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T 
BR  =    ■=- 

1 


(2) 


Therefore  it  is  clear  that  if  the  ratio  of  ten  sion  to  current  in  the  wire 
is  the  same  as  the  ratio  of  momentum  to  charge  of  the  particle,    the 
wire  will  assume  the  path  of  the  charged  particle. 

The  above  derivation  is  for  a  constant     B,    but  it  is  possible  to 
divide  the  wire  into  many  small  segments  where  the  immediate  field 
is  considered  constant,    and,    in  the  limit,    the  wire  orbit  and  particle 
path  will  be  the  same  for 


T_    =    P 

l  e 


(3) 


In  practice  the  wire  should  be  very  fine  and  flexible  so  that  the 
force  of  gravity  may  be  neglected.      In  workable  units, 


...       /    .         2.94  M(g) 
P(Mev/c)=       i(amp)  ' 


where     M    indicates  the  weight  attached  to  the  wire  to  provide  the 
tension. 


2,      Example 

Tie 
point  /         Wires 


Pulley 


Fig.    16.      Experimental  setup  in  wire -orbit  method. 
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A  typxcal  experimental  setup  for  using  the  wire -orbit  method 
is  shown  in  Fig.    16.      The  two  pieces  of  wire  passing  through  the 
magnet  are  in  series  so  that  they  represent  two  possible  "rays"  of 
particles  eminating  from  the  tie  point,    A.     Actually  the  tie  point  must 
be  located  slightly  ahead  of  the  planned  point  of  particle  emission, 
since  the  wires  occupy  physical  space  and  cannot  be  made  to  emanate 
from  a  point.      The  wires  then  pass  through  the  magnet,    past  a  frame 
for  measuring  their  separation,    and  then  to  a  point  just  short  of  a 
pulley.     A  short  piece  of  insulated  wire  is  tied  to  the  current-carrying 
wire  and  then  passed  around  the  pulley.      This  is  done  so  that  the 
currents  in  each  wire  have  the  same  direction  sense.      Using  the  tie 
point  as  the  object  and  the  intersection  of  the  wires,    as  determined 
by  a  projection  of  the  measured  separation  at  the  frame  and  diameter 
of  the  pulley  (as  in  Fig.    17)  as  the  image,    one  solves  the  focusing 
problem  by  simple  geometry.     In  practice  the  current  in  the  magnets 
is  adjusted  until  the  separation,    s,    at  the  frame  is  correct  to  provide 
the  desired  image  point. 


s  =  separation  at  frame 
d  =  pulley  diameter 


Fig.    17.      Geometric  determination  of  focus. 
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mation, apparatus,  method,  or  process  disclosed  in 
this  report. 


As  used  in  the  above,  "person  acting  on  behalf  of  the 
Commission"  includes  any  employee  or  contractor  of  the  Com- 
mission, or  employee  of  such  contractor,  to  the  extent  that 
such  employee  or  contractor  of  the  Commission,  or  employee 
of  such  contractor  prepares,  disseminates,  or  provides  access 
to,  any  information  pursuant  to  his  employment  or  contract 
with  the  Commission,  or  his  employment  with  such  contractor. 


